QFT results for neutrino oscillations and New Physics 



D. DelepineEl and Vannia Gonzalez Macia^ 
Physics department, Division de Ciencias e Ingenieras, Universidad de Guanajuato, 
Campus Leon, CP. 37150, Leon, Guanajuato, Mexico. 

Shaaban KhaliHI 

Centre for Theoretical Physics, The British University in Egypt, 
El Sherouk City, Postal No, 11837, P.O. Box 43, Egypt 

G. Lopez CastrcH 

Departamento de Fisica, Cinvestav, Apartado Postal 14-740, 07000 Mexico D.F., Mexico 

(Dated: April 30, 2009) 

The CP asymmetry in neutrino oscillations, assuming new physics at production and/or detection 
processes, is analyzed. We compute this CP asymmetry using the standard quantum field theory 
within a general new physics scenario that may generate new sources of CP and fiavor violation. Well 
known results for the CP asymmetry are reproduced in the case of V -A operators, and additional 
contributions from new physics operators are derived. We apply this formalism to SUSY extensions 
of the Standard Model where the contributions from new operators could produce a CP asymmetry 
observable in the next generation of neutrino experiments. 
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I. INTRODUCTION. 

Since the experimental results implying that the neu- 
trinos are massive Q , it is expected to have CP violation 
phases in the leptonic sector. However, the situation in 
the leptonic sector is very different from the quark sector 
where the CP violation is clearly established in K and B 
mesons physics. The only evidence for flavor violation 
in the leptonic sector comes from neutrino oscillations 
and there is, so far, no confirmation for CP violation in 
leptonic decays. Hence, measuring any CP asymmetry 
in neutrino oscillation will open a new window to study 
CP violation and related problems as leptogenesis or CP 
violation in r decays. New physics beyond the Stan- 
dard Model (SM), like low energy supersymmetry, may 
be probed at the LHC and consequently new sources of 
CP and flavor violation can be expected. These new 
sources of CP and lepton flavor violation, also classified 
as non-standard interactions (NSI), could give important 
contributions to the CP asymmetry in neutrino oscilla- 
tion (see ref.@,illi]). 

In addition, a proper procedure to take into account 
these non-standard interactions is important as neutrino 
experiments are reaching a high level of accuracy. Even 
if neutrino masses and lepton mixing matrices solve the 
solar and atmospheric neutrinos anomalies, it is clear 
that NSI could affect the oscillations parameters as de- 
termined by the next generation of neutrino experiments. 
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Effects of NSI have been widely studied, considering their 
contributions to solar and atmospheric neutrino prob- 
lems, in neutrinos factories, in conventional and upgraded 
neutrino beta beams, e^e~ colliders, neutrino-electron 
and neutrino- nucleus scattering and in many other as- 
pects of neutrino physics 0, 0, H, i, i, [M [Hf E, [11 • 

The interpretation of a possible CP asymmetry in neu- 
trino oscillation due to the SM or NSI is still an open 
question. The following two hypothesis are usually con- 
sidered in the analysis of CP asymmetry in neutrino os- 
cillations: (i) The probability of a process associated to 
neutrino oscillation can be factorized into three indepen- 
dent parts: the production process, the oscillation proba- 
bility and the detection cross section, (ii) The CP asym- 
metry in this process is due to the CP violating phase 
in the lepton mixing matrix. In a pioneering work, the 
authors of Ref. [1] have studied CP violating effects due 
to contributions from new neutrino interactions in the 
production and/or detection processes in neutrino oscil- 
lation experiments. However, only corrections to the V-A 
SM charged current interactions were considered In 
Ref. 3, the (V-A) (V-A) and (V-A)(V-hA) operators as- 
sociated to muon decays, but not to the pion decay, have 
been considered. In muon decays, the interference be- 
tween SM and New Physics of (V-A)(V-|-A) operators 
is suppressed by |30i |. In Ref. |I4{ . the quantum 

field theory formalism has been used to describe neutrino 
oscillations but New Flavour Interactions were not taken 
into account. 

The goal of this paper is to go beyond this approx- 
imation and to propose a generic framework based on 
quantum field theory to get a simple expression for the 
CP asymmetry without imposing any assumptions on the 
operators generated by New Physics. We shall show that 
in such a case, new contributions to the CP asymme- 
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FIG. 1: Feynman diagram of the process 7r(pi) — > jJ.'^{p2) + 
//^(p) followed by the detection process:i';(p) + N{pn) 
N'{pj^i) -\- l{pi). It is important to recall that the flavour 
of neutrino intermediary state is unobservable. 

try appear and it could be important to take them into 
account once we want to constraint new physics using 
experimental data. We shall illustrate this in the case of 
the supersymmetric extension of the Standard Model but 
it is clear that our treatment is valid beyond that and it 
can be applied to NSI effects in all neutrino experiments. 

II. GENERAL FORMALISM 

Let us start by giving a sketch of the idea of the present 
work. Let us consider a virtual neutrino that is produced 
at the space-time location {x,t), travels to (x , i ) and is 
detected there because it interacts with a target produc- 
ing a charged lepton /. For definiteness, we illustrate this 
process with the production of the neutrino in tt"*" decay 
and its later detection via its weak interaction with a 
target nucleon N (see Figure 1): 

^ yf{p) + N{pn) ^ iV'bjv') + Kpi) ■ 

We shall call i'*'^ state, respectively the neutrino which 
is produced at source in conjunction of a /i^ and the 
neutrino which is detected through the detection of a 
charged lepton of flavour I. These effective states are not 
necessary of or I flavour once NSI are introduced. 

Energy-momentum conservation at the production and 
detection vertices requires pi = P2 + P and p + pn = 
Pn' + Pi- Weak interactions in the SM acting at the 
production and detection vertices conserve lepton num- 
ber and this process is interpreted as a flavor change: 
a /i-neutrino is transformed into a /-neutrino due to os- 
cillation. However, in the presence of new physics flavor 
violating weak interactions can occur, and the flavor iden- 
tification of the neutrino at the production and/or decay 
locations via its associated charged lepton is not unique 
anymore. 

This process is relevant for neutrino factories and long- 
baseline accelerator experiments such as superbeams. Pi- 



ous are produced at the source when a high-energy pro- 
ton beam (~ 10^ GeV) hits a target (made of solid dense 
material e.g. Be, Al, Ca, etc.). Superbeam neutrino 
sources are mainly from pion decays while neutrino fac- 
tories sources are mainly from muons. 

If we assume that light neutrinos are left-handed, Lep- 
ton Flavor (LF)-violating semi-leptonic interactions can 
be described by the following effective Hamiltonian: 

n = 2V2GFVad{C^{jirPL'^k){ujaPLd) 

+C\ {-pPLi^k) (uPRd) 

+C^ii i-pao^fsPLUk) {ua'^^PRd) 

+C^L{P^c.pPLVk){uCJ^''PLd)] , (1) 

where k runs over the three leptonic flavors and Pr l — 
(l±75)/2. 

In the following and for simplicity, we consider the 
case where LF violation occurs only at the 7r+ decay ver- 
tex. It is straighforward to include the effects of such 
New Physics at the detection vertex using this formal- 
ism. Note that the tensor currents proportional to the 
^5L{R) Wilson coefficients will not contribute to 7r+ decay 
because it is not possible to generate an antisymmetric 
tensor from the pion momentum alone. Thus, the only 
non-vanishing hadronic matrix elements at the produc- 
tion vertex are: 

(0|d7^75uk+> = ^^p'i (2) 

(0M75«k+) = , (3) 

mu + md 

where /^r = 130 MeV is the pion decay constant and ^ 
denote the light quark masses. 

Using the relation i^k — Uka^a between flavor k and 
mass a neutrino eigenstates, we get the following ampli- 
tude for 7r+ decay: 

{^i+vk\n\'^+)^u{p)o''v{p2) , (4) 

where: 

V2 \mu + md 

-i{Cl-C',)K) (5) 

with = j^pl^. 

Now consider a neutrino of flavor k that is produced 
in 7r+ decay via LF-violating interactions and is detected 
with flavor I at a later time via (LF-conserving) charged 
current scattering off the nucleon N (Fig. 1). In Quan- 
tum Field Formalism, the S-matrix amplitude for the evo- 
lution for the system from initial state is given by: 

= |d4^dV5:e^(P'-^"+^«')-'^(J^^,)M 
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where A'-^ {x — x) basically describes oscillation of neu- 
trinos during its propagation. It is important to stress 
that we never introduce Energy-momentum eigenstates 
to describe the neutrinos as in our formalism the neu- 
trinos appear as virtual particles. In terms of massive 
neutrino propagators we can write 



(2^ 



-ip{x —x) 



nil,. + ie 



where Un are the usual elements of the Umns niixing 
matrix. We write this propagator in a more convenient 
form by integrating upon the time component of the four- 
momentum [ISl j: 



Ai\x - X) =Y,UuUt, j ^ 

-iE{t —t)^i~p{x 



(6) 



{E^^.j" + m^je{t -t) 



(-i;^..7° - -p.^ + m^,)e{t - t) 



where E^. — ^Jp^ + m^. . As usual, we interpret the first 
term as neutrinos propagating forwards in time and the 
second as anti-neutrinos propagating backwards in time. 

Thus, by keeping only the first term of the propagator 
in Eq. (6) we get the amplitude 



GpVud 



75 J 



(7) 



with T — [t — i) > is the time elapsed from the pro- 
duction to the detection space-time locations of neutri- 
nos and Pp = pi — Pn + Pn' ■ Equivalently, the time- 
dependent amplitude from initial to final states is the 
integrand of eq.(7): 

ns-.^{T) = {27r)^S^iPF+P2-pi)iGFVud)HJNN')^ 



x{m^Al + B*^^)v{p2)e'P-^ 



where 



2E,, 



Ak 



niu + nid 



(8) 



(9) 
(10) 



The time evolution amplitude for the corresponding 
CP-conjugate process which correspond to the observa- 
tion at source of a fi^ and detection of a is given by: 

rp._p,(T) = {27r)H\PF+P2-pi){GFVudf{JNj,')^, 
'Xf7r^Vi{PlW\^-l5)ih~ i>N+ i)N') 



X [m^Ak + Bk^^) u{p2)e'^''^ 



2E,, 



(11) 



In order to get the CP asymmetry in the general case, 
let us choose the following form of the nucleon weak ver- 
tex: 

(JNN')t^ = Ujv'(pjv')7m(5v +5A75)uw(P7v) (12) 

with gv = gviq^ = 0) = 1 and gA = gA{q'^ = 0) « -1.27 
[iq . Under these approximations, one has: 

= ^{B*m^ - m^A*){Bkm^ - mT,Ak) 

q.k 

C//*C/fe,[/zjC/^e^"(^'^"-^''.)F(P, M) ,(13) 

where F{P, M) is a kinematical function that depends 
on masses and momenta of external particles but not on 
neutrino flavour and will drop in the acp{T) asymmetry 
defined as: 

, , |r.g-.(r)p-|r,;,-.,(r)P 



|T.j_,(r)P-MT^._^,(r)P 
D{r) 



(15) 



III. CP ASYMMETRY FROM Ci WILSON 
COEFFICIENT ((V-A)(V-A) OPERATOR) 



In the usual formalism developed by refs. 0]0 one 
considers only New Physics corrections due to = 
CsAii^kfj. + <^fik)- In this approximation, one gets for the 
amplitude 

r,._,,(T) = {2n)*e'P°^^S\PF+P2-pi){GFVud)'{JNM'),^ 

x/7rC'5MUi(pi)7^(l - 75)(A" ^>N+ ^>n') ^>-^v{P2) 



2E^^ 



(16) 



Once taking the |Tya_y,(T)p, the first two lines of the 
previous equation will give us the kinematical part of the 
process which is common to all neutrino flavor and the 
last line will give, at first order in neutrino masses, the 
usual neutrino oscillation results including Non Standard 
Interactions. It is interesting to note the equivalence of 
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our expression with the approach done in ref.[3| where 
they defined: 

K E (f^eVpe + C/;,(l + ep^) + U;,e^r) > (17) 



where > are the neutrino mass eigenstates and {v^ > 
is the initial flavour state produced at neutrino source. 
Computing one obtains the neutrino oscil- 

lation part of |Ti/s_i^j once assuming that \/2Ey. — 
{\/2E„){l + 0{ml.)) and keeping the first term of the 
expansion in rn^. . 



IV. CP ASYMMETRY FROM Csa WILSON 
COEFFICIENTS 

The interesting results of this formalism is that natu- 
rally all NSI can be taken into account without making 
any a priori assumptions. In particular, one can easily in- 
clude the effects of new CP violating phases and flavour- 
violating interactions from Cs^i Wilson coefficients which 
involve scalar and pseudoscalar density operators. As an 
example, let us assume that the only sources of CP vio- 
lating phases are coming from the scalar operators C^^ 
(equivalently a relative weak phase between and Bi^, 
different from rnr), so the numerator of the CP asymme- 
try will be proportional to 



N{t) oc i?e , ^ 




AlB.UuUlMjUqji 



-iT{E„ 



The Standard model contributes only to ~ CsMSkfi+ 
corrections from New Physics (in our convention, Csm — 
1) and C3.4 are produced by New Physics. At first order 
in new physics, one has to replace Bq by (5,^ in the CP 
asymmetry as previously defined. 

It is interesting to note that we could get this result 
by assuming from the beginning that 



coefficients C3 4 as mentionned previously in refs.[T3| [Sl- 
it is clear that the limit on helicity of the muon in pion 
decays [11] [III will constraint the contribution of (73,4 not 
to be bigger than a few percent. It is important to note 
that most of the constraints on e.y come from four Icpton 
Fermi Operators (see ref.[l|) as in — > eee or rare muon 
decays as /i — > ej; also some lepton flavour violating tau 
decays impose very strong constraints on e^i with i = e 
or T. In the case of operators which contribute to pion 
decays, one should note that these operators can be writ- 
ten as the product of a leptonic current and a hadronic 
current. Thus, the Wilson coefficents in both cases can 
not be related to each other unless they can be factor- 
ized into the product of a leptonic and a hadronic part 
and if CP and flavor violations occur only in the lep- 
tonic vertices. This procedure is model-dependent and 
should be verified in each specific model of New Flavour 
Interactions. 



V. 



APPLICATION TO SUSY MODELS 



Let us apply this formalism to the usual Minimal Su- 
persymmetric extension of the SM. Hence, as an appli- 
cation, we calculate the supersymmetric (SUSY) contri- 
butions to the dominant pion decay mode tt^ fJ-^Vf^ 
(or its charge conjugate) [13, [H, [2^ through the Wilson 
coefficients and focus our attention on the SUSY constri- 
butions to scalars and pseudoscalars Wilson coefficients 
(C3 and C4). 



The effective lagrangian contributing to the i^^ 
flavor change at the neutrino source, is given by 



acp{T) 



\{^iKi^W~\mK{r)) 



\{^i\K{r)W + m\K{r))\ 



(18) 



where one deflnes the 



with 



E 



(f^eV^e + U;,{1 + e^^) + U;,€^r) W > (19) 



Csm 



It is worth mentioning that using the 



QFT formalism, all e^j^ are expressed in terms of the Wil- 
son coefficients and once the e^f are defined, the analy- 
sis considered in Ref.[l,|l| can be easily implemented. As 
one can see from eq.(ll), the New Physics is enhanced by 
a factor mT^/imu + md) « 15 which multiplies the Wilson 



where Cj are the dimensionless Wilson coefficients and 
Oj are the relevant local operators at low energy scale as 
defined in Eq.(2). The leading order contributions to C3 
and C4 under the experimental constraints for the Wil- 
son coefficients induced by SUSY are explicitly shown in 
the expressions below. These describe box type diagrams 
of chargino-neutralino exchanges. Other SUSY contribu- 
tions (vertex corrections) are suppressed by the Yukawa 
couplings of light leptons. In fig ([2]), the Feynman dia- 
grams of dominant SUSY contribution to the pion decay, 
7r~ /i~ + represented. 

To simplify the expression, we shall assume that 
the squark and slepton masses are degenerated {fhf ~ 
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FIG. 2: Feynman diagrams of dominant SUSY contribution 
to 7r~ /i~ + iv^. 
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MNS 



XM-^O^M-^-J{Xa,Xb) 



X,Y,A,B 



ZVZVud^F -J 



^^^^^^^ t/A2C/fy^*M 



7-7 ttMNS 
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where the N and U, V matrices are the matrices which re- 
spectively diagonaUze the neutraUnos and chargino mass 
matrices, g is the gauge weak couphng and 9w is the 
Weinberg angle. The sum on X, Y, A, B is respectively 
on the squarks, sleptons, charginos and neutralinos. The 
fimctions /(xa, a;;,), J{xa, a^b) are given by 



I{Xa,Xb) 



J [Xa 1 Xb) 



1 



1 



IGir^nif^ \xa — Xb, 



1 



1 



IGn'^nif'^ \Xa ~ Xb, 

- 1 -j- Xfi Xa In Xa 



and Xa b are defined 



{Xa Xb) 



(Xa Xb) 



In or- 



der to have non-vanishing asymmetry we should have a 



non- vanishing imaginary part of C3 or C4. The com- 
plex phase of C^^^ can be due to the matrices that di- 
agonalize the squark mass matrix or the chargino mass 
matrix. If we assume that the squark matrices are di- 
agonal, i.e the diagonalizing matrices are identity, we do 
not have any source of CP violation from squark ma- 
trices. Then the only remaining source is the phase of 
fi which induces a complex phase in chargino mass ma- 
trix and hence U and V unitary matrices. This phase 
is strongly constrained by the neutron electric dipole 
moment (EDM) [HI, [H, ^ [13, [H, [H but it is pos- 
sible to avoid this constraint if the first generation of 
squarks are heavy enough and decouple from the two 
other squark generations. Another way to avoid the neu- 
tron EDM constraint is to assume flavor violating struc- 
ture in the squarks and sleptons mass matrices. This 
allows to have new sources of CP violation which could 
contribute to any CP violating observables and/or could 
relax the constraint on the fi term coming from neutron 
EDM limit [23]. 

Assuming that the CP violating sources comes from 
the fj, term in the chargino sector, e^e is given by C^^^ 
times the enhancement factor mj;/ [rriu + md)- For typ- 
ical chargino and neutralino masses of order 150 GeV 
and sfermion masses around 100 GeV, it is possible to 
obtain e^e as large as 10"*^. In order to maximize the 
absolute values of e's, one assumes that charginos and 
(20) neutralinos have quasi-degenerated masses. In figure ([3]), 
we present our numerical results for |e^e| (solid line) 
and |e^T-| (dashed and dotted lines) as function of the 
chargino mass, for tan/3 = 50 and tan/3 = 10. As ex- 
pected from the expressions of e^e is not sensitive 
to tan /3 if 6f^^^ ~ 0. One should emphasize that the 
SUSY model presented in this section is the simplest one 
and we could expect enhancement in flavour violation 
effects once a non-universal structure is assumed in the 
soft-SUSY breaking terms. In this case, the values of 
or/and e^r may be significantly enhanced. A detailed 
analysis for the CP asymmetry in neutrino oscillation in 
SUSY model with non-minimual flavor will be considered 
elsewhere. 

As one can see from Fig. ([3]), Cg^ is typically of order 
lO"'^. These values lead to CP asymmetry in neutrino 
oscillation of order 10~^ ~ 10~^, depending of the pa- 
rameters of the neutrino beam and the size of the exper- 
iment baseline Such asymmetries are reachable at 
next generation of reactor and beam neutrino oscillation 
experiments 0, [1, [3, 1, [1] ■ 



VI. CONCLUSIONS 

To conclude, we proposed a formalism based on Quan- 
tum Field Theory where it is possible to include all 
sources of Non Standard Interactions including new CP- 
and flavour-violating interactions. We show that using 
this method, it is straightforward to include the effects 
of the scalar and pseudoscalar operators densities which 
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FIG. 3: Absolute value of e^e (solid line) and e^jr (dashed and 
dotted lines) for different choice of tan /3, respectively tan /3 = 
50 and 10. The other parameters are given hy fhf — 100 GeV 
and m^o = 110 GeV. 



corrections to (V-A)(V-A) Wilson coefficients. Tliis ap- 
proach was well justified when neutrino experiments were 
not so accurate and other sources of New Physics could 
be easily neglected. But now neutrino experiments enter 
in the field of high precision experiments where the effects 
of NSI coming from (S-P)(S+P) and (S-P)(S-P) opera- 
tors could have an impact on the determination of lepton 
mixing parameters including the CP violating phases. To 
illustrate our results, we applied it to SUSY models where 
we assume that the only CP-violating phases appear in 
^3,4 Wilson coefficients which correspond respectively to 
(S-P)(S-P) and (S-P)(S-|-P) operators. In such a case, we 
show that with reasonable values for the SUSY parame- 
ters, it is possible to generate a CP asymmetry as large 
as 10-1 - 10-^ 



appears in any New Physics Models. In the limit where 
we consider only (V-A)(V-A) operator, we reproduce 
the usual results reported in Ref.jl]. It is important 
to emphasize that most of the studies on NSI in neu- 
trino physics have been done assuming that New Physics 
contributions to neutrino interactions is mainly due to 
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